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ABSTRACT 


The  free-electron  model  for  conjvigated  systems  is  consistently  developed  as  the 
limiting  case  of  a three-dimensional  q\iantvim-mechanlcal  treatment  of  the  7r-electrona 
in  such  systems.  Joint  conditions  (for  branching  points)  and  boxmdary  conditions  (for 
free  endpoints)  are  derived,  and  the  hermit icity  of  the  Hamiltonian  is  shown.  A matrix 
formulation  of  the  theory  is  established  which  makes  the  application  to  large  systems 
feasible  and  at  the  same  time  leads  to  a close  analogy  with  the  LCAO  model  (>  LCAO  MO 
treatment  considering  only  nearest  neighbor  interactions).  Quantities  analogous  to  the 
quantities  q (=  charge  in  an  atomic  orbital)  and  p (»  bond  order)  are  defined,  and 
special  attention  1s  given  to  alternant  conjugated  systems,  for  which  a population 
theorem,  analogous  to  the  one  in  LCAO  theory,  is  valid. 


INTRODUCTION 


SYSTEMS  of  conjugated  double  bonds  in  organic  molecules  are,  in  a good  approximation, 
described  by  assuming  that  the  6-electrons  form  bonds  which  maintain  the  molecular 
frame  and  that  the  7r-electrons  move  in  the  potential  of  this  frame.  Interesting  pro- 
perties of  the  molecules  depend,  in  a first  approximation,  on  the  w-electrons,  e.g. , 
the  relative  bond  lengths,  certain  types  of  chemical  reactivity,  ultraviolet  si>ectra, 
and  (to  a lesser  degree)  dipole  moments;  and  it  is  therefore  of  interest  to  investigate 
the  notion  of  the  Ti-electrons  by  itself.  A basic  characteristic  of  this  motion  is  the 
high  mobility  or  delocalization,  i.e . , the  fact  that  the  different  electrons  are  not 
bound  to  definite  bonds,  but  move  rather  freely  over  the  whole  conjugated  structure. 

The  description  by  molecular  orbitals  (MO's)  seems  therefore  particularly  suited  for 
this  case. 

The  classical  method  of  constructing  MO's  consists  of  forming  them  as  LCAO  MO's 
(linear  combination  of  atomic  orbitals).  The  rigorous  application  of  this  method  leads 
to  an  extremely  large  nurber  of  difficult  integrals,  and  only  benzene  has  been 

^Thls  work  was  assisted  in  part  by  the  Office  of  Naval  Research  under  Task  Order  IX 
of  Contract  N6orl-20  with  The  University  of  Chicago. 
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BuooeaafuLly  treated  In  this  way.^  In  order  to  achieve  a more  general  method  appli- 
cable to  larger  claaaes  of  molecalea,  one  generally  makes  the  asBtmptlon  that  only 
nearest  neighbors  interact  and  oonsldera  the  remaining  Integrals  as  parameters  in  the 
final  result.  This  simplified  method^  haa  been  rather  successful,  and  in  the  following 
we  always  mean  this  simplified  LCAO  theory  whenever  we  mention  the  LCAO  model  without 
further  specif loatlon. 

On  the  other  hand,  the  existence  of  the  core  potential  suggests  the  possibility  of 
expressing  this  potential  in  terms  of  coordinates  which  are  fitted  to  the  particular 
frame  of  the  molecule,  e.g^. , cylindrical  coordinates  in  the  case  of  benzene.  As  a 
consequence,  the  molecular  orbitals  would  be  expressed  in  terms  of  the  same  ooerdl- 
nates.  Again,  this  approach  can  be  carried  through  for  benzene,^  but  it  becomes  very 
complicated  for  larger  moleoules.  Here  too,  however,  there  exists  a simplified  formu- 
lation applicable  to  any  conjugated  systemt  the  one-dimensional  free-eleotron  model, 
This  model  is  not  new,^  but  it  has  mostly  been  regarded  as  a plotorial,  but  rough, 
description.  It  is  the  purpose  of  the  present  paper  to  show  that  the  free-eleotron 
model  deserves  the  same  recognition  and  standing  as  the  simplified  LCAO  model.  A close 
analogy  between  the  mathematics  of  the  two  models  will  be  derived,^®  and  in  the  second 
paper  of  the  present  series^  a close  agreement  will  be  established  between  the 

^M.  Ooeppert-Mayer  and  A.  L;  Sklar,  J.  Chem.  Phys.  6,  645  (1958))  C.  C.  J.  Roothaan, 
and  R.  Q.  Parr,  J.  Chem.  Phys . 17,  1001  (1949) i Parr,  Craig,  and  Ross-,  J.  Chem.  Phys. 
18,  1561  (I950)j  0.  C.  J.  Roothaan,  Rev.  Mod.  Phys.  2J,  69  (1951). 

®See,  e.g. , Coulson  and  Longuet -Higgins,  Prop.  Roy.  Soo.  A191.  59  (1947))  see  also 
reference  8. 

^See  K.  Ruedenberg  and  R.  0.  Parr,  J.  Chem.  Phys.  Ig,  1268  (1951). 

Pauling,  J.  Chem.  Phys.  4,  675  (1936))  0.  Schmidt,  Ber.  Deutsoh.  Chem.  Qes.  73A, 

97  (1950);  J.  R.  Platt,  J.  Chem.  Phys.  1^,  484  (1949);  H.  Kuhn,  Helv.  Chim.  Acta  31. 
1441,  1780  (1948),  2247  (1949),  24,  1508,  2571  (1951),  J.  Chem.  Phys.  16,  840 

(1948),  II,  1198  (1949),  Zeit.  f.  El.  Chem.  §2,  I65  (1949),  2^,  220  (195D) 

S.  Nlkltlue,  J.  Ch^.  Phys.  6l4  (1950). 

^*’n.  S.  Bayliss,  J.  Chem.  Phys.  I6,  287  (1948),  II,  1855  (1949),  Aust.  J.  Scl.  Res. 

109  (1950);  N.  S.  Bayliss  and  J.  C.  Riviere,  Aust.  J.  Sol.  Res.  A4,  344  (1951)) 

W.  T.  Simpson,  J.  Chem.  Phys . 16,  1124  (1948). 

^°An  analogy  between  the  LCAO  model  and  the  PE  model  in  the  simplest  cases  (no  branch 
points)  wan  already  noticed  by  W.  T.  Simpson,  J.  Chem.  Phys.  II,  1218  (1949),  and 
H.  H.  Jaffe,  J.  Chem.  Phys.  20,  1646  (1952). 

5c.  W.  Scherr,  J.  Chem.  Phys.  21,  000  (1953). 
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numerical  resulte  obtained  by  the  two  modela.  Thus,  from  an  empit'loal  point  of  view, 
the  free-eleotron  model  seems  to  serve  as  well  to  desoribe  conjugated  eyatems  as  the 
simplified  LCAO  model. 

But  also  from  the  theoretical  point  of  view,  the  free-eleotron  model  offers  sever- 
al attractive  properties;  l)  once  the  free-eleotron  model  is  chosen,  the  quantum- 
mechanical  problem  is  rigorously  solved)  one  makes  no  further  approximations  of  virtu- 
ally uncontrolled  Implications,  as,  e.£.,  the  assumption  of  considering  only  nearest- 
neighbor  interactions  in  the  LCAO  model)  S)  in  the  case  of  hydrocarbons,  the  free- 
eleotron  model  has  only  one  parameter,  the  neighbor  distance  D,  as  compared  to  the  two 
parameters  B,  S of  the  LCAO  model )^  3)  the  meaning  of  the  parameter  D in  terms  of 
measurable  quantities  is  more  precisely  defined  than  the  meaning  of  the  parameters 
occurring  in  the  LCAO  model)  therefore  there  la  but  little  “leeway"  to  "adjust"  D and 
the  free-eleotron  model  comes  closer  to  being  an  "absolute  theory"  (a  no  adjustable 
parameters))  4)  the  free-electron  wavefunotlons  are  easy  to  vlsuallze--a  fact  which  is 
demonstrated  by  the  model  of  J.  R.  Platt  explained  in  the  third  paper  of  this  series 
Of  course  the  LCAO  model  has  the  advantage  that  it  is  in  a position  to  desoribe  simul- 
taneously all  electrons,  not  only  the  n-eleotrons,  in  a molecule. 

The  second  part  of  the  present  paper  is  restricted  In  scope  to  the  case  that  all 
atoms  in  a conjugated  structure  are  mathematically  equivalent,  l..e. , practically  to  the 
case  of  hydrocarbons.  A method  to  Include  the  effect  of  heteroatoms  is  in  preparation. 

1 . GENERAL  CONCEPTS 

a.  The  Model 

The  organic  molecules  to  which  the  present  theory  applies  are  characterized  by 
the  fact  that  each  of  them  contains  a so-called  conjugated  system  of  bonds,  jL.e. , a 
system  of  alternating  single  and  double  bonds.  Each  atom  within  a conjugated  system 
contributes  two  fl-electrons  and  one  w-electron  for  the  chemical  binding  of  the  conju- 
gated structure.  In  the  approximation  considered  here,  the  6-electrons  form  6-bonds 
which  maintain  the  geometrical  arrangement;  and  the  molecular  framework,  stripped  of 
the  7T-electrons,  forms  a "skeleton"  or  "core"  which,  because  of  positive  charge, 

^Por  a fair  comparison  It  should  be  noted  that  the  correlation  of  the  LCAO  model  with 
experimental  data  can  be  obtained  by  assigning  a fixed  value  to  5 and  treating  p as 
the  only  empirical  parameter.  The  LCAO  model  actually  contains  a third  parameter,  o, 
which  however  serves  only  to  calibrate  the  zero  point  of  the  energy. 

R.  Platt,  J.  Chem.  Phys . 2^,  000  (1953). 
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creates  a potential  under  whose  Influence 
all  TT-electrons  move  throughout  the  entire 
conjugated  structure.  Thus,  If  formula  I 
of  Figure  1 describes  for  example  the 
naphthalene  molecule,  formula  II  In  the 
same  figure  wculd  be  an  appropriate  de- 
scription of  the  corresponding  molecular 
skeleton.  Indeed,  this  formula  gives  even 
a correct  geometrical  picture  of  the  physi- 
cal skeleton,  since  It  Is  well  known  that 
all  atoms  belonging  to  a conjugated  system 
lie  in  one  plane.  , 

The  starting  point  of  the  theory  to  be  developed  here  Is  the  assumption  that  the 
core  potential  Is  infinitely  high  everywhere  except  on  the  bond  lines  of  the  skeleton 
structure  II,  where  It  Is  finite.  Obviously  suoh  a potential  has  to  be  considered  as 
the  limiting  case  of  a potential  which  Is  infinitely  high  everywhere  except  within  a 
thin  tube  extending  along  the  b'^nd  lines.  Furthermore  the  Intereleotronlc  Interaction 
is,  In  a first  approximation,  neglected  between  the  Tr-eleotrons  so  that  the  latter 
move,  independent  of  each  other  under  the  influence  of  the  core  potential, 
b . Linear  Bond  Paths 

In  order  to  determine  the  Implications  of  the  assumptions  Just  made,  let  us,  as  a 
first  example,  consider  the  motion  of  an  electron  in  a long  thin  box  of  length  I and 
with  a square  cross-seotlon  of  side-length  e « i (see  Figure  2)  This  model  corre- 
sponds to  the  TT-electronio  motion  In  a polyene,  e.g,. , butadiene.  The  eigenfunctions 

*nmm' = (2/l)^(2/e  )sin(iTnx/l)sln(7imy/e)sin(7rm'a/e)  (1.1) 


Fig.  1 
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must  bp  anti0ymmet;’lc  with  reapeot  to  the  plane  z a e/2,  which  we  take  to  be  the  nodal 
plane  of  the  n-eleotrons  ao  that  the  quantum  numbera  asaume  the  values 

n ® 1,2, ?,••■}  m ■ 1,2,7, •••}  m'  " 2,4,6, ■■■  (1,2) 

From  the  exproaaion  for  the  eigenvalues 

- (hVSm)  [(n/l)^(mA)®+(m»A)®l  , (1.3) 

It  follows  that,  because  of  e « £,  the  motions  in  the  y-  and  z-dlreotiona  remain  in 
their  ground  state  (m  b i,  m<  <>  2)  while  the  x-motlon  ia  excited  to  ita  lower  levels, 
so  that 

E - a (hV8m)[  (n/^)®+(5/e®))  “ E^+  conat.  , (1.4) 

E^  - (h®/8mi®)n®  , (I.5) 

and  the  wavefunotions  have  a single  nodal  plane  perpendicular  to  the  z-axla  in  agree- 
ment with  the  nodal  properties  of  Bprr  atomic  orbitals.  Hence  the  y-  and  z-motlon  may 
be  neglected  for  low  exoitationa,  and  it  is  sufficient  to  work  with  the  energies  (1.5) 
and  the  wavefunotions 

0n(x)  “ Bin(7rnx//)  , (1.6) 

which  may  be  defined  by  the  equations 

<t>n(x)  » / dx  i dy4>®^g(x,y,z)  = ^ 

sign  of  <fj^(x)  B sign  of  «i„jg(x,yBie,ZBie)  , (1.7') 

where  Q denotes  a cross-section  of  the  tube  and  dQ  = dxdy,  a cross-section  element. 

The  functions  (1.6)  are  the  normalized  eigenfunctions  of  the  one-dlmenaional  eigenvalue 
problem  considered  by  the  free-eleotron  model,  and  the  energies  E^  (1.5)  are  the  corre- 
sponding energy  eigenvalues. 

The  limiting  process  e -*  0 leaves  0jj(x)  and  unaffected  and  therefore  meaning- 
ful, whereas  and  become  Infinite.  The  latter  circumstance  means  that  the  y- 

and  z-excltation  energies  become  Infinitely  high  so  that  and  give  (for  6 = 0) 

a correct  description  of  the  system  for  all  finite  excitations  of  the  x-motlon,  !•£. , 
for  all  finite  quantum  numbers  n. 

Since  high  excitations  effect,  however,  structural  changes  in  organic  molecules. 
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they  do  not  have  to  he  oonaldered  and  the  aaeumption  e « / la  actually  aufflclent  to 
Justify  the  of  the  one-dlmenalonal  deaorlption  given  by  (1.5,6),  l.e.,  the  physical 
assumption  la  only  that  the  Tr-eleetrona  form  a oloud  whose  length  la  large  compared  to 
its  width.  The  passage  to  the  limit  e -♦  0 la  only  a matter  of  mathematical  oonvenlenoe 
and  does  not  Imply  that  the  Tr-elootron  oloud  la  aotually  Infinitely  thin, 
c . Qomparlaon  with  the  LCAO  Model 

The  present  theory  is  thus,  physically,  a three-dimensional  theory,  and  the  elec- 
tron cloud  of  a TT-eleotron  la  Ita  MO.  These  MO's  are  not  entirely  different  in  nature 
from  the  MO's  formed  as  linear  combinations  of  atomic  orbitals  (LCAO  MO's),  if  all 

matrix  element a between  non-neighbor Ina  atoms  are  neglected- -the  assumption  adopted  In 
o 

the  LCAO  MO  model.”  Indeed,  such  an  assumption  implies  that  the  radii  of  the  (finite) 
atomic  orbitals  are  smaller  than  one  bond  length  In  straight  chains  and  emaller  than 
iVj  n 0.87  bond  lengths  in  benzene  rings,  so  that  the  molecular  orbitals  formed  by  them 
obtain  also  a tubular  shape  as  shown  in  Figure  3,  where  the  circles  indicate  the 
atomic  orbitals,  and  the  shaded  area  is  inaocesslble  to  the  Tr-eleotrons, 

The  MO's  constructed  as  LCAO  MO's  correspond  to  a potential  which  Is  not  constant 
within  the  accessible  region.  This  generalization  can  be  readily  Introduced  In  our 
model  by  assuming  a potential  V ° V(x)  inside  the  box  of  Figure  2.  Thereby  the  eigen- 
values (1.5)  hnd  the  eigenfunctions  (1.6)  will  be  changed,  but  they  will  remain  finite 
so  that  all  arguments  presented  above  remain  unaltered. 

There  Is  no  need  to  Introduce  a potential  whioh  depends  also  on  y and  z,  since, 
because  of  e « t,  the  motions  in  the  y-  and  z-dlrections  will  in  any  case  stay  In  the 
ground  state  so  that  the  (y,z) -dependence  remains  irrelevant.  The  LCAO  MO  model  (with 
neglect  of  non-neighbor  interaction)  also  considers  essentially  only  "longitudinal 
excitations",  since  one  can  see  from  Figure  3 that  the  "transversal  profile",  l.e.,  the 
(y,z) -dependence,  of  the  LCAO  HO  la  much  less  changed  by  varying  the  coefficients  of 
the  different  atomic  orbitals  than  the  "longitudinal  profile". 

The  purpose  of  the  foregoing  remarks  is  not  to  Justify  the  free-electron  model  by 
vague  analogies  with  the  LCAO  MO  model,  but  to  give  a plausible  reason  for  the  agree- 
ment which  exists,  as  we  hope  to  show,  between  the  two  methods. 


®An  Introduction  to  the  LCAO  MO  treatment  and  further  literature  may  be  found  in  B.  and 
A.  Pullmanns  Les  theories  Electron iques  de  la  chlmle  organlque  (Masson  and  Co., 

Paris:  1952),  ch.  IV,  2. 
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d.  Non-Linear  Bond  Paths 

Let  us  takie  benzene  as  an  example  of  the  case  that  the  bond  path  Is  no  longer  a 

straight  line,  and  let  us  for  the  moment  assume  that  the  Tr-eleotrons  move  In  a thin 

olroular  tube  (torus  of  radius  R)  rather  than  In  a hexagonal  tube,  Let  us  further 
replace  the  annular  well  potential  representing  the  torus  by  the  analytic  potential 

V(r,d,z)  - U(2a)‘®{(Z/R)2+[(r/R).(R/r)]2}  , (1.8) 

U - fiVsiPR^  > (1.9) 

where  the  cylindrical  cocrdlnates  z,  r,  and  0 (z  perpendicular  to  the  molecular  plane) 
are  used,  and  a Is  a dimension  leas  constant  « 1.  It  Is  sensible  to  define 

€ =.  aR  (1.10) 

as  the  diameter  of  the  corresponding  tube,  whereas  & = 2ttR  la  Its  length,  so  that 
Indeed  e/£  = a/27r  « 1 . 

It  has  been  shown^  that  the  potential  (1.8)  leads  to  energy  levels 


n_ 
0 r z 


(1.11) 


where  E corresponds  to  angular  excitations,  E corresponds  to  radial  excitations, 

"r 

and  E corresponds  to  out-of -plane  excitations.  The  energy  levels  E of  the  "longl- 
'^2  "0 
tudlnal  notion"  are  given  by 
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= Un^f(a^n|)  , n.  - 0,  *1,  ±2,  ...  , (1.12) 

f(x)  o ((l+4x)i-l]/2x  » 1 - X + ...  , (1.13) 

and  for  the  two  transversal  motions  one  finds  the  excitation  energies 

dE„  ° AE„  - 2U/a  » . (l.l4) 

“r  "g  "(j 

It  Is  therefore  again  true  that  the  excitations  of  the  transversal  motions  can  be 
neglected,  and  that  the  energy  levels  of  the  longitudinal  motion  approximate  the 
limiting  levels 

.Urn  E„  - Un?  . (1.15) 

6-0  ^ 

These  are,  however,  Identical  with  the  levels  obtained  by  Considering  the  one- 
dlmenslonal  motion  of  a particle  along  the  line  which  Is  the  limit  of  the  tube  as  e— 0. 

These  conclusions  can  be  expected  to  be  Independent  of  the  particular  example 
chosen,  l.e.,  to  hold  not  only  for  a circular  tube,  but  for  a tube  along  any  arbitrary 
bond  path,  e.g. , for  the  hexagonal  tube  of  benaene.  This  expectation  will  be  Justified 
at  the  end  of  the  next  subsection.  We  have  then  the  general  result  that  the  assumption 
of  a potential.  Infinite  everywhere  except  on  the  bond  skeleton,  is  equivalent  to 
describing  the  electronic  motion  by  an  one-dlmenslonal  MO,  <^(x),  depending  on  the 
"bond-path  coordinate"  x,  and  assuming  for  ^(x)  the  Schriidlnger  equation 

{(d/dx)2+(2m/ft®)[E-V(x))}0(x)  = 0 ^ (l.l6) 

where  V(x)  Is  an  appropriate  potential  along  the  bond  skeleton,  The  complete  descrip- 
tion of  the  electron  la  then  furnished  by  a three-dimensional  molecular  orbital  which 
can  approximately  be  written  as 

*(x,y,z)  « d>(x)(2/e)sln(7ty/€)sln(27rz/6)  , (1.17) 

0 s y S:  £ , -|e  « z « |€  , £ « £ . (1.17') 

Here  the  x-ooordlnate  la  tangential  to  the  bond  paths,  y Is  perpendicular  to  x In  the 
molecular  plane,  and  z Is  perpendicular  to  the  molecular  plane,  which  Is  chosen  to  be 
the  plane  z * 0.  The  orbital  4>(x,y,z)  Is  antisymmetric  with  respect  to  this  plane.  A 
more  rigorous  form  for  <t(x,y,z)  la  obtained  by  solving  the  thln-tube  problem  corres- 
ponding to  the  bond  skeleton;  between  such  a <J>(x,y,z)  and  the  one-dlmenslonal  orbital 
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0{x)  onfi  has,  .in  analogy  with  (l.7»7')»  the  relation 

0^(x)  « 11m  f dQ4>^(x,y,z)  , 

€-0  Q(x) 

sign  of  0(x)  = sign  of  <^(x,y•€/2,z•»e/4 ) , 


(1.18) 

(1.19) 


where  e la  the  diameter  of  the  tube,  and  Q(x),  its  oross-seotlon  at  the  point  x. 
e.  Branching  Points  (Joints) 

Those  points  In  a conjugated  structure  where  three  branches  meet  (we  call  them 
branching  points  or  Joints)  present  the  one -dimensional  formulation  with  a particular 
problem.  Let  us  split  the  skeleton  of  Figure  1 (II)  up  Into  three  branches  and  define 
on  each  branch  an  Independent  arbitrary  coordinate  system  as  shown  by  Figure  1 (III). 
Let  ^^(xg)  be  the  part  of  the  wavefunotlon  d>(x)  which  lies  on  the  branch  B(nl,2,5)j 
more  precisely)  let  *g(xg)  be  defined  by: 


<?b(xb)  ='  0(x),  If  X ■ Xg,  l.e.,  on  branch  B, 
0b(Xb)  ® 0,  If  X # Xjj,  l.e,,  on  ail  other  brarohes, 


(1.20) 


so  that 


0(x)  - E 0B(Xg)  . (1.21) 

B 

Then  the  Schrfidlnger  equation  (1.16)  holds  for  each  branch  function  0g(xg),  and  the 
question  arises!  how  to  fit  together  those  (tgCx^)  which  meet  In  one  Joint. 

A first  condition  Is  quite  clearly  the  "continuity  condition":  the  three  branch 

functions  0j^(xj^),  0g(x2),  0j(x^)  must  assume  the  same  value  at  their  meeting  point, 
l.e ! 


*^1^*1J^  “ ' (1.22) 

where  Xgj  Is  the  coordinate  of  the  Joint  on  branch  B. 

However,  since  Eq . (I.l6)  is  a second  order  differential  equation,  we  have  to 
derive  also  a condition  involving  the  first  derivatives  of  the  functions  Od(Xb).  For 
this  purpose  we  consider  again  the  system  of  thin  tubes  corresponding  to  a bond  skele- 
ton. Let  Figure  4a  Illustrate  the  part  of  this  system  in  the  neighborhood  of  a Joint. 

As  before,  let  4>(x,y,z)  be  the  solution  of  the  three-dimensional  Schrildinger 
equation  inside  the  tube,  fulfilling  the  boundary  condition  4>(x,y,z)  « 0 at  the  walls. 
From 
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(b) 

Pig.  4 

hence  It  follows  that  also 


div(<{>  grad  0)  « grad  0 + , 


(1.25) 


we  obtain  by  virtue  of  Oauss ' theorem 


dV  (grad®*  + *£i*)  => 

where  S la  the  surface  of  the  volume,  V,  and  (3/Sn)  la 
the  derivative  in  the  direction  of  the  outer  normal  on  S. 
For  the  volume,  V,  we  choose  the  volume  Indicated  by  the 
shaded  area  in  Figure  4b.*  Since  ♦ vanishes  on  the  walls, 
and  Qj,  Qg,  are  planar  cross-sections,' Eq . (1.24) 
simplifies  to 

I ^ dQ<J>®  = 2J  , (1.25) 

H 

J = ^ dV(grad®4>  + 4>A4>)  , (1.25') 


I dS*(S*/&nj  = dS(a*®/an)  (1.24) 


whore  the  coordinate  X|^  Is  tangential  to  the  tube  branch 
B,  and  pointing  away  from  the  Joint.  Now  consider  first 
the  limit  e -•  0 (c  = tube  diameter).  By  virtue  of  Eq, 
(1.18, 20)  the  left-hand  side  of  Eq.  (1.25)  goes  Into 


£ 

B 


’‘b“^BJ+6 


= finite;  (1.26) 


11m  J = finite, 

G-»0 

and  therefore  the  Integrand  of  J becomes  infinite  of  the  order  e'®  as  e -<  0.  Secondly, 
consider  the  limit:  {llm  e -•  0 and  11m  6 -»0},  where  6 Is  the  distance  of  each  of  the 

three  surfaces  Qg  from  the  Joint  midpoint  (see  Figure  4b).  Since  the  surfaces  Q.^  are 
purely  fictitious,  and  * does  not  fulfill  any  boundary  conditions  on  them,  the  limit 
6 -.0  (in  distinction  from  the  limit  e -*  0)  does  In  no  way  cause  * to  become  Infinite. 
Hence,  In  the  limit  {e  -« 0,  fi  -*0},  the  Integrand  of  J,  once  again,  goes  to  infinity 
as  e'®,  and  hence 


llm  J = 0 . 

4 -O,  e—O 


Therefore,  the  expression  (I.26)  must  vanish  for  Xg  Xgj,  and  we  find  the  new 
boundary  condition 
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f " I ° *B  “ *BJ  '. 

which  by  virtue  of  the  continuity  condition  (1.22)  Is  equivalent  to 


(1.S7) 


^B  Joint  “ ® 

Since  Eq.  (1.28)  may  be  derived  also  in  the  case  of  complex  wave functions,  it  is 
clear  that  Eq.  (1.28)  expresses  the  conservation  of  the  quantum-mechanical  current 
densltyi 


(lfi/2m)[0(d«*/<lx)  - ^*(d0/dx)]  . (1.29) 

However,  since  all  wavefunetions  0(x)  which  occur  as  solutions  of  our  problems  are 
(with  a few  exceptions)  real,  their  current  density  vanishes  everywhere  ao  that  the 
current  conservation  becomes  trivial.  We  shall  call  Eq.  (1.28)  the  "conservation 
oondltlon." 

The  derivation  of  the  conservation  conditions  is  Independent  of  the  number  of 
branches  meeting  at  one  Joint.  Although  for  our  purposes  only  the  case  of  three 
branches  per  Joint  Is  useful,  it  is  still  of  interest  to  note  that  also  in  the  general 
case  the  conditions  (1.22)  and  (1.28)  suffice  precisely’ to  detemine  the  solution  of 
the  eigenvalue  problem,  if  the  wavefunotlon  is  required  to  have  a unique  value  every- 
where. In  order  to  recognize  this  independence,  consider  a system  oontaining 

a^  Joints  oolleotlng  ? branches 
a^^  Joints  collecting  4 branches 
• « * 

aj  Joints  collecting  J branches, 

and  b "free  endpoints." 

If  B___  is  the  total  number  of  branches,  then 
max 

2Bmax  = 3aj  + 4aj^  + ...  + Jaj  + b . (l.30) 

Since  Eq.  (l.l6)  Is  a second-order  differential  equation,  the  left-hand  side  of  Eq. 
(1.30)  gives  the  number  of  constants  In  the  total  wavefunotlon  0(x)j  the  right-hand 
side  of  Eq.  (1.30),  on  the  other  hand,  gives  the  total  number  of  Joint  and  endpoint 
conditions.  Since  there  is  one  more  parameter  (the  eigenvalue)  and  one  more  condition 
(the  normalization  of  the  wavefunotlon)  in  the  problem,  the  number  of  unknowns  to  be 
determined  equals  the  number  of  equations  to  be  fulfilled. 
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It  IB  furthermore  easily  seen  that,  in  the  case  of  one  closed  branch,  as  In 
benzene,  the  Joint  conditions  are  equivalent  to  the  condition  that  the  wavefunotlon 
have  a period  that  Is  an  integral  fraction  of  the  branch  length. 

On  the  basin  of  the  Joint  oondltlons  (1.22,28)  It  Is  possible  to  give  a more 
rigorous  Justification  for  the  statements  leading  up  to  Eq.  (l.l6),  namely,  that  there 
la  no  difference  In  the  treatment  of  linear  and  non-linear  bond  paths.  Consider  any 
non-branching  point  on  a bond  path  as  a Joint  with  two  branches;  the  conditions  (1.22, 
28)  then  are  applicable  and  yield  the  result  that  both  0(x)  and  (d^/dx)  are  continuous 
at  any  such  point.  These  continuities,  however,  entail  the  validity  of  the  Schrfidlnger 
equation  (1.16)  throughout  the  entire  bond  path  containing  only  such  points.  Consider 
in  particular  a "corner"  in  the  bond  path,  Si£.,  in  naphthalene  [see  Figure  1 (II)]:  on 
the  two  straight  portions  of  the  bond  path  meeting  at  the  corner,  Eq.  (I.I6)  le  valid; 
at  the  corner,  0 and  (d0/dx)  have  to  be  continuous;  and.  If  V(x)  is  continuous,  Eq. 

(1.16)  yields  that  the  left  limit  of  (d/dx)^0  equals  the  right  limit  of  (d/dx)®0.  It 
follows  that  (d/dx)  0 exists  at  the  corner  and  equals  the  common  limit  so  that  Eq. 

, (1,16)  Is  valid  at  the  corner.  Hence,  corners  in  a bond  path  have  no  implications,  If 
the  TT-eleotronlc  potential  1s  continuous. 

The  conditions  (1.22,28)  were  first  used  by  Hans  Kuhn.^  As  yet  no  proof  or  dis- 
cussion of  the  conservation  condition  (1.28)  seems  to  have  been  given, 
f . Formulation  of  the  Oeneral  Mathematical  Problem 

It  may  be  useful  to  collect  at  this  point  the  equations  which  determine  the  one- 
dlmenslonal  free-eleotron  molecular  orbital  (PEMO).  They  are: 

1)  the  Schrddlnger  equation  on  the  branches: 

H 0(x)  = {(-fi^/2m)(d/dx)^  + V(x))0(x)  = E 0(x)  , (1.31) 

2)  the  Joint  conditions: 

[0j^(xj^)  » •<>g(x2)  “ at  Joint  (continuity  condition)  , (1.32) 

£ (S0g/5Xg) = 0 (conservation  condition)  , (1.33) 

B 

3)  the  boundary  conditions  for  free  endpoints i 

0(x)  = 0 at  free  endpoint.  (1-34) 

Kuhn,  Helv.  Chlm.  Acta  32.  224?  (1949). 
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Coneernlng  the  branches  having  free  endpoints,  there  Is  still  a question  that  remains 
to  be  aettled--namely , where  to  locate  the  endpoint.  For  reasons  which  will  become 
clear  later  [see  appendix]  we  adopt  the  following  postulate:  The  free-electron  paths 

terminate  one  bond  length  beyond  the  last  atom.  Thus,  for  the  styrene  molecule  of 
Figure  5a,  whose  bond  skeleton  is  given  by  Figure  5b,  the  free-electron  path  Is  pic- 
tured In  Figure  5ot  it  extends  up  to  the  point  E,  whereas  the  last  atom  Is  situated  at 

The  definition  Initially  given  for 
the  free-electron  model  has  therefore  to 
be  modified.  It  should  now  read:  The 

potential  for  the  7r-eleotronlo  motion  Is 
assumed  to  be  infinite  everywhere  except 
on  the  tree-electron  path,  where  it  is 
finite  j the  free-ele,dtrOn  path  Is  Iden- 
tical with  the  bond  skeleton  except  that 
It  exceeds  the  latter  at  its  free  end- 
points by  one  bond  length.  Some  authors 
have  used  different  conventions  concerning 
the  free  endpoints  In  order  to  fit  certain 
experimental  dataj  our  choice  is  based  on  theoretical  grounds,  as  will  be  shown  in  the 
append lx . 

£,.  Hermltlolty.  Orthogonality.  Normalization 

The  problem  just  formulated  Is  an  eigenvalue  problem  in  a one-dlmenslonal,  but 
multlconnected  configuration  space.  Let  [see  Eqs.  (1.20,21)] 

f(x)  « E fsCxg)  (1.35) 

B 

be  a function  In  this  spaoe  satisfying  the  Joint  and  endpoint  conditions  (1.32,33,34), 
but  otherwise  arbitrary.  The  Integral  over  the  configuration  apace  Is  then  defined  as 

/ dxf(x)  = E / dXj,f(x„)  , (1.36) 

& B ig  “ 

where  I indicates  the  total  length  of  the  free-electron  paths  and  ig  Indicates  the 
length  of  branch  number  B.  In  connection  with  such  integrals  the  following  generaliza- 
tion of  Green's  theorem  is  important. 

Consider  the  integration  over  the  one-dlmenslonal  space  Indicated  in  Figure  6, 


the  point  C. 


(a) 


(b) 

Pig.  5 


(C) 


8 
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where  the  coordinates  Xg  runs  from  0 (at  the  Joint)  to  ag  (at  the  endpoints).  Let 
f(x)  and  g(x)  be  two  functions  of  the  type  (1.35)  which,  however,  do  not  vanish  at  the 
points  Sg,  since  Figure  8 Is  assumed  to  be  part  of  a larger  bond  skeleton.  Then  the 
following  theorem  holds : 

Indeed,  Green's  theorem  yields 

^ dXg[fg(d/dXg)®gg  - gg(d/dXg)^fg]  = (fg(d/dXg)gg  - gg(d/dXg)fg]^  (1.38) 

for  eaoh  branch.  If  one  now  sums  over  the  three  branches,  those  terms  on  the  right- 
hand  Bide  of  (1.38)  which  are  taken  at  the  Joint  (Xg  a o)  vanish  by  virtue  of  the  Joint 
eondltlons  (1.32,33,34) i 

S [fB(d/dXg)gg  - SQ(d/<lXg)fg]j|,  Q “ f(0)  E (dgg/dXg)g  - g(0)  E (dfg/dXg)Q  = 0 
B B B B 

and  hence  Eq.  (1.37)  results. 


Pig.  6 


If  one  now  decomposes  the  total  free-electron  path  of  a bond  skeleton  into  parts 
of  the  form  of  Figure  6 and  applies  (1.37)  to  eaoh  part,  then  all  right-hand  side 
contributions  of  Eq.  (1.37)  cancel  each  other  or  vanish  at  the  free  endpoints,  where 
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/ dxf(x)(d/dx)^g(x)  H / dxg(x)(d/dx)^f(x)  , (l.?9) 

I i 

where  the  Integration  extends  over  the  total  bond  slceletoni  i.e.,  (d/dx)‘^  Is  a her- 
mltian  operator  with  respect  to  our  multloonnected  configuration  space,  and  so  Is  the 
Hamiltonian  H of  the  Sohrfldlnger  equation  (I.31).  As  a consequence,  the  eigenvalues 
of  that  equation  are  real . and  eigenfunctions  belonging  to  different  eigenvalues  are 
mutually  orthogonal.  (Eigenfunctions  Belonging  to  the  Bume  BigBnValus  can,  of  course, 
be  chosen  so  as  to  be  mutually  orthogonal.)  The  validity  of  these  properties  depends 
essentially  upon  both  joint  conditions  (1.52,33). 

The  eigenfunctions  are  required  to  be  normalized! 

/ dx  = r / dx-0i(xo)  •>  1 . (1.40) 

t B 

2.  CASE  OP  CONSTANT  POTENTIAL 
.£.•  Oeneral  Form  of  the  Solution 

If  the  potential  V(x)  Is  constant  throughout  the  entire  free-electron  path,  It 
can  bo  assumed  to  vanish  identically,  and  the  SchrBdinger  equation  (I.31)  Is  readily 
solved.  The  eigenfunction ‘of  the  energy  level  Ej^  has  the  form 

0^(x)  = S * (2*1) 

s 

“ *^n°®®(‘*nV*Bn^  ' (2.2) 

with 

E^  = (fiV2m)k2  . (2.3) 

In  the  following  we  shall  omit  the  subscript  n wherever  It  Is  not  needed.  The  condi- 
tions which  determine  the  parameters  k,ag,6g  (b  = 1,  2,  3,...)  become  very  simple.  The 
continuity  conditions  (I.32)! 

a^co3(kXjj+6^)  = agCOs(kXgj+62)  = ar^oos  (kXjj+6^)  ; (2.4) 

the  conservation  conditions  (1.33) s 

a^sin(kXjj+6^ ) + a2Sln(kXgj+6g)  + ajSln(kXjj+6j ) « 0 , (2.5) 


or,  by  virtue  of  (2.4): 
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tan(kx^j+6^)  + tan(kXgj+6g)  + tan(kx^j+6j)  « 0 . (2.6) 

In  these  equations,  x^j  denotes  the  coordinate  of  the  Joint  on  the  branch  B.  The 
normalization  condition  becomes 

i S = / dx  <l>^(x)  = 1 , (2.7) 

B I 

Where  ig  denotes,  as  us\ial,  the  length  of  the  B-th  branch.  Eq.  (2.7)  is  established  in 
the  following  way.  On  branch  B one  obtains 

g 2 2 ^B 

^ dXg  a|oo8  (kXg+6g)  = i-ag  ^ dXgtl+oos  2(kXg+6g)]  , 

= ia|{£g  + k"^[sln(kX3+6g)eos(kXg+6g)]Q®}  , 

jL  I 6 e j 

ig 

^ dXB03(xB)  “ ia|^3  + ik”2{[<ji3(d^)3/dn) ]q  + («g(di^g/dn) ] , (2.8) 

where  (d/dn)  denotes  the  derivative  at  a Joint  (or  endpoint)  In  that  direction  on  the 
branch  which  points  away  from  the  Joint  (or  endpoint).  By  virtue  of  the  normalization, 
the  left  hand  side  of  (2.8)  will  become  unity  after  summation  over  all  branches}  hence 

k2{2  - Saf^g}  = S {[^B(d«g/dn)]Q  + ^ ^ 

B B B 

Here  each  branch  enters  on  the  right-hand  side  only  by  its  first  and  last  point.  These 
are  either  Joints  or  free  endpoints.  The  latter  do  not  contribute  anything  since  0 >=  0 
at  a free  endpoint.  The  terms  arising  from  Joints,  on  the  other  hand,  may  be  ordered 
according  to  Joints,  so  that  the  right-hand  side  of  (2.6)  may  be  written 

S ‘bJi  joint  J ■ 

where  B = 1,  2,  3 denote  the  three  branches  merging  in  Joint  J and  E denotes  the  sum 

J 

over  all  Joints  of  the  bond  skeleton.  By  virtue  of  Eqs . (1.32,33),  each  term  in  the 
summation  over  J vanishes  separately,  and  thus  Eq.  (2.7)  Is  proved.  This  equation  may 
be  considered  a generalization  of  the  familiar  result  obtained  in  the  case  of  a single 
branch,  namely  that  the  mean  value  of  cos  t between  0 and  tt  Is 

It  Is  useful  to  Introduce,  at  this  point,  a few  conventions.  First  of  all  we  note 

that 


33 


RUEDENBERQ  AND  SCHERR 


agCos(kXg+6g)  - a'gOOB(kXg+6'g)  , 

with 

a'e  “ -ag  , 6'g  «■  6g  ± 7T  . 

In  order  to  eliminate  any  ambiguity  we  shall  always  ohoose  the  phase  suoh  that 

iw  > fig  > ' (2.11) 

thereby  simultaneously  fixing  the  relative  signs  of  the  amplitudes  Sg  of  the  different 
branches.  At  the  same  time  this  establishes  a one-to-one  oorrespondence  between  6g  and 
tanfig. 

There  Is,  of  course,  still  a factor  e^"*'  (7  ■»  real)  which  la  arbitrary  in  the 
normalized  functions  (2.1).  .However,  we  require  the  wavefunotlon  to  be  real,  and  thus 
only  a factor  (-1)  remains  arbitrary  in  the  total  function  d(x)  on  Eq.  (2.1). 

In  the  second  place,  it  la  seen  that  tne  sign  of  k and  6 may  be  simultaneously 
reversed  without  affecting  the  wavefunotlon  (2.2)  or  the  eigenvalue  (2. 5).  We  adopt 
the  convention  always  to  ohoose 

k 0 (2.12) 

which  entails  In  each  case  a definite  sign  for  6.  The  de  Broglie  wavelength  of  the 
motion  described  by  the  molecular  orbital  (2.1)  is  then  given  by 

- 27r/k  . (2.13) 

.‘Jlnce  a constant  potential. V(x)  treats  all  atoms  as  equal,  it  is  natural  to 
consider  the  distance  D between  neighboring  atoms  as  a constant  throughout  the  conju- 
gated structure.  We  adopt  this  assumption  and  define  the  dimensionless  constant 

K =■  IcD  , (2.14) 

which  will  prove  to  be  convenient.  Eq.  (2.3)  can  then  be  written 

(En/En)  = (a/D)2|c2  ^ (2.14') 

where  a = ■ft^/e^m  Is  the  Bohr  radius,  and  Eg  = 6^/23  is  the  ionization  potential  of  the 
hydrogen  atom.  It  will  soon  become  clear  that  for  our  purposes,  k lies  always  in  the 
range 
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0 « K « 7T  , (2.15) 

so  that  a one-to-one  oorreapondenoe  exlata  between  /c,  sins,  and  cosk. 

Elgenveotora  of  a Free-Electron  Eigenfunction 
Let 


^(1),  41(2),  ...  «{N)  , (2.16) 

denote  the  valuea  which  the  eigenfunction  'Pix) , (2<1),  aaaumes  at  the  N atoms  of  the 
conjugated  structure.  We  define  the  N-dlmensional  row  vector 


0(1) 

0(8) 


I 0(N)  J 


(2.17) 


and  call  It  the  free-eleotron  eigenvector  corresponding  to  the  eigenfunction  0(x), 
(2,1),  Thlfl  vector  la  uniquely  determined  by  the  function  0(x);  however,  Inversely 
0(x)  Is  also  uniquely  determined  by  the  eigenvector  0 In  conjunction  with  Its  eigen- 
value (2,J).  Indeed,  each  branch  B contains  at  least  two  atoms,  say  P and  Q,  so  that 
the  constants  a^  and  can,  by  virtue  of  (2.2),  be  expressed  In  terms  of  0(P),  0(Q), 
and  k ■ (2mE/ft^)^.  In  order  to  obtain  explicit  expressions  for  a^  and  6g  it  is  con- 
venient to  employ  the  relation 

i0g(xg+5)  = oos(k§)0g(xg)  + k"^8in(kO(d0g/dXg)jj  , (8. 18) 

B 


where  | Is  arbitrary.  Eq.  (2.18),  which  will  also  prove  to  be  useful  later,  follows 
from  (2.2)  by  virtue  of  the  addition  theorem  for  the  cosine  or  else  by  means  of  the 
Taylor  expansion 

0B(xg+§)  ■=  £^(?Vn! )d”0g(Xg)  = e^^0g(Xg)  , 

(>b(Xb+5)  = (cosl4d  - 1 slnU3)0B(xg)  , (2.19) 


where  3 = (d/dx).  The  equivalence  of  (2.18)  and  (2.19)  follows  from  the  Schradlnger 
equation 

(d/dx)^(?.(x)  + k^$(x)  =.  0 , (2.20) 


In  practice  one  chooses  the  origin  of  the  coordinate  on  a branch,  always  such  as 
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to  have  It  coincide  with  an  atom  or  with  the  midpoint  of  a bond.  Consider  therefore 
first  the  case  that  Xg  = 0 coincides  with  an  atom  P so  that 

^g(P)  - Sgcoafig  , (2.ai) 

and  let  Q denote  the  atom  next  to  P at  Xg  = D,  so  that  by  virtue  of  (2.18)  and  (2.14) 


<ig(Q)  =>  $g(P)[oosK  - Blnit  tanfig] 
From  (2.21)  and  (2.22)  one  obtains  the  relations 


(2.22) 


tandg  - [oosK  - 0g(<l)/03(P)I/8lnic  , (2.23) 

s|  - 0|(P)/coB®6g  B (1  + tan®6g)(^3(P)  , (2.24) 

O 

for  the  calculation  of  tan&g  and  Sg.  The  sign  of  Sg  must  be  the  same  as  the  sign  of 
0g(P),  since  in  Eq.  (2.21)  one  has  oosflg  > 0 because  of  the  convention  (2.11),  6n  the 
other  hand,  If  Xg  *»  0 coincides  with  the  midpoint  of  the  bond  between  the  atoms  P and 
Q,  then  application  of  (2.18)  for  Xg  =>  0,  4p  «•  -iD,  and  yields 

liiB(P)  ■»  ag[ooB^K  eos6g  + Bin^K  slnfig]  , (2.25) 

^b(Q)  ■ ag[oosiK  oosdg  - slnix  sinfig]  , (2.26) 

or 


♦g(P)  + 0g(Q)  = EagCosix  cos6g  , 
dg(P)  - $b(Q)  *>  EagSln^K  sinfig  . 

P 

Hence  one  obtains  for  tanfig  and  agS 

1 + oosx  $b(P)  - 0g(Q) 

" " sThx  - 5g(P)'  '+  *g(Ci)  ' 

a|  «■  (0|(p)  + 0|(Q)  • 2008K  0g(P)0g(Q)]/8ln®K 


(2.27) 

(2.28) 

(2.29) 

(2.30) 


since  oosfig  > 0,  the  sign  of  ag  is  the  same  as  the  sign  of  [0g(P)  + 0g(Q)],  see  Eq. 
(2.27). 

It  Is  useful  to  define  another  N-dlmenslonal  row  vector,  namely 


^ 8 


01 

’2 

ft* 

0M 


(2.31) 


r- 
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whose  elements  tpg,  . . . have  subscripts  corresponding  again  to  the  N atoms  of  the 

system  and  are  defined  In  the  following  manner  (see  2.16) s 


9p  " D'^d(P)  , If  P la  not  a Joint  ("ordinary  atoms")  , 
B (D5/2)^<i(Q)  , If  ft  la  a Joint  . 


(2.52) 


The  vector  f has,  as  we  shall  show,  the  property 

" A ’p  ” 5 J 

P-1  '"D  i Big  " 

Where  denotes  the  hermltlan  conjugate  matrix  of  which  here  is  identical  with  the 
transposed  matrix.  Since  ^(x)  is  assumed  to  be  a normalized  eigenfunction  (see  1.51), 
we  obtain 


* N 
^ ^ « S 

P-1 


(2.34) 


and  therefore  call  f the  normalized  free  electron  eigenvector. 

Zn  order  to  establish  (2.55)  consider  first  the  n atoms  ft^,  ftp,  ft^,  ...  on  one 
branch  B.  If  Xp  Is  the  coordinate  of  the  first  atom  ft^^  • Pg  and  Xg  ■ Xp  + (n-l)D  Is 
the  coordinate  of  the  last  atom  ■ Eg,  then  let  us  denote  the  argument  of  the  cosines 
in  (2.2)  by  (see  2.14) 


f - kXp  + flg  (first  atom  Pg) 

e - kXg  + 6g  - f + (n-l)x  (last  atom  Eg)  . 


(2.55) 


For  the  sum 

N p 

3-  - » ♦l(Qv)  » (2-56) 

v-1 

one  obtains 

Sb  - a|  cos^(KV+f)  » ia^  £ {2  + e^2(*v+f)  g-12(icv+f)j 

= ia^(2n  + e^2f g-12f j 


and  by  virtue  of  (2.35) 


Sg  - ia|(2n  + - e"^^®"^‘'’''h/[e^'‘‘-e"^''] ) 

Sg  - i^a^in  + ^[sln(2e+K)  - sin(2f-«) J/clnx} 
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Sb  “ ia|{(  n-l)  + ootc(8lne  cose  - alnf  ooaf)  + ooa^e  + ooa^f}  , (2.37) 

Qlvlng  ig  thB  same  meaning  aa  In  Eq.  (2.7)  and  giving  (d/dn)  the  same  meaning  aa  In 
Eq.  (2.8)  one  oan,  In  the  case  of  a branch  between  two  Joints,  write  (2.37)  as 


DSg  - ia|/g  + iD[«|(Pg)  + + (D/2k)ootK{[$3(d0g/dn)]p^  + [«<ig(d((»g/dn)  ]j^} 


(2.38) 


If  a branch  B contains  a free  endpoint  and  is  the  last  atom,  then  th*  free  end- 
point with  the  boundary  condition  0 - 0 Is  located  one  bondlength  beyond  e . , at 

the  site  of  a fictitious  atom  Since  ^(0^+1^  " 

n o o 

Sg  (n)  = d|(ft^)  - d|(Q^)  - SB(n+l)  j 

and  we  obtain  Eq.  (2.38)  agalnj  but  all  terms  coming  from  the  free  end  vanish. 

Let  us  now  sum  Eq.  (2,38)  over  all  branches  of  a conjugated  structure  and  write 
the  result  In  the  form 

D fi  {Sg  - id§(Pg)  - ^|(Eb)} 

■ X ia|ig  + (D/2k)oot*:  S {[«B(d0B/dn) ]p^  + [0g(d(<iB/dn) ]j^}  . (2.39) 

faking  Into  consideration  that  all  contributions  from  free  endpoints  vanish,  one 
recognizes  that  by  virtue  of  (2.36)  and  (2.32)  the  left-hand  side  of  (2.39)  Is 
Identical  with 


No# 

P-1  ^ ^ 


where  the  sum  over  P Includes  all  atoms  of  the  conjugated  system.  Furthermore,  the 
second  term  on  the  right-hand  side  of  (2.39)  is  of  the  same  nature  as  the  right-hand 
side  of  (2.9),  and  vanished  for  the  reasons  given  there.  Thus  Eq.  (2.39)  becomes 
identical  with  Eq.  (2.33),  and  the  proof  Is  completed. 

The  normalization  theorem  is  not  limited  to  the  eigenvector  f (2. 31).  Choose,  on 
the  free  electron  path,  any  set  of  L-equldlstant  points  Rj,  Rg,  ...  Rj^,  whose  neighbor 
distance,  d.  Is  adjusted  In  such  a manner  that  all  Joints  and  endpoints  are  contained 
in  the  set,  and  define  the  L-dlmenslonal  row  vector 


M 

f - 


£1 

’2 


(2.40) 


I 
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~ 1 

a d^0(R)  , If  R is  not  a Joints 

B (d3/2)^0(R)  . if  R is  a Joint. 


(2.41) 


Prom  the  preceding  derivation  It  Is  apparent  that  the  normalization  thnorem  (2.42,43) 
holds  equally  for  ^ . 

For  example,  one  may  ehooae  d ■ ^D,  and  then  the  points  Rj,  Rg,  R^  are  all 
atoms  Plus  all  bond  midpoints.  One  obtains 


~2  i_2 

“ *’r 


at  atoms i R » P , 


hence 


5^  ■ iD0‘^(R)  at  bond  midpoints*  R = M , 


r 5p  • i r ipp  + iD  z 0®(H) 

R ” P M 


(2.42) 


By  virtue  of  the  normalization  theorem,  the  left-hand  side  of  (2.42)  is  unity,  and  the 
first  term  on  the  right-hand  aide  la  hence  one  finds 


D Z 0‘'(M)  B 1 , 
K 


(2.43) 


where  the  sum  extends  over  all  bond  midpoints.  Eq.  (2.43),  which  may  be  called  the 
normalization  theorem  for  the  bond  midpoint,  serves  as  a convenient  numerical  check  in 
any  computation  involving  the  0®(H). 

A different  example  is  provided  by  the  stilbene  molecule  (see  Figure  7).  There 
it  is  possible  to  choose  a four-dimensional  normalized  row  vector  corresponding  to 
the  four  atome  A,  B,  0,  D,  with  the  neighbor  distance  d = 3D. 


Fig.  7 
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0 . Matrix  Form  of  the  Eigenvalue  Problem 

Consider  three  ornseoutlve  points  P_j,  Pq,  on  a branch  with  coordinates 

(Xg-D),  Xg,  (Xg+D)  respectively.  Applying  the  relation  (2.18)  for  0g(p_j^)  and  0g(P^j)» 
one  obtains 


«g(P_^)  - 200BK  ^g(PQ)  + ^b(P+i>  “ 0 . (2.44) 

If  Pq  Is  the  last  atom  at  a free  end,  then  ♦(P.j)  =•  0,  and  hence 

-2C0BK  ^(Pq)  + «b(P+i)  “ 0 (2.45) 

In  that  ease.  Consider  furthermore  a Joint  J and  Its  three  neighboring  atoms  P2,  Pg, 
Pj.  Applying  Eq.  (2.18)  for  (P^)!  (Pg)»  (Pj)»  and  adding  all  three  relations  one 
obtains 

3 3 

E *o(Pb)  - 300SK  d(J)  + k"^Bln/e  E j • 

Bal  Bal 

Because  of  the  conservation  condition  (1.55) > the  second  term  on  the  right-hand  side 
vanlshoa,  and  the  result  Is 

0(Pl)  + d(Pg)  + «(pj)  - 500BK  0,(J)  n 0 . (2.46) 

Slnoe  for  each  atom  one  has  one  relation  of  the  type  (2.44,45,46),  one  finds 
exactly  as  many  equations  as  there  are  atoms.  Hence  one  has  a system  of  N homogeneous 
linear  equations  for  N quantities  , 0(Pg),  ...  0(Pjj).  This  may  be  written 

P 0 « 0 , (2.4T) 

where  0 Is  the  eigenvector  defined  In  Eq.  (2.17),  and  P Is  a matrix  whose  structure  la 
most  easily  explained  by  writing  It  down  In  a special  case,  say  styrene.  If  the  atoms 
are  numbered  In  the  way  Indicated  In  the  free-eleotron  path  of  Figure  5c,  then  the 
matrix  assumes  the  forms 
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where  all  elements  which  are  omitted  vanish,  and 

F H 2cosk  . 

Let  T be  the  diagonal  matrix  defined  by 


'^PQ  “ *PQ  ’ 

1,  if  P is  not  a Joint  , 

’ 1 2 4 

(j)*j  if  P is  a Joint  , 

so  that  the  normalized  eigenvector  (2. 31, ^2)  is  given  by 

^ a D^T"^«  . 


(2.50) 


(2.51) 


Then  Eq.  (2.47)  oan  be  rewritten  as  the  eigenvalue  equation 

F#  - P#  , (2.52) 

where  F la  the  matrix 

F - TUT  + PI  , (2.53) 

(X  a unit  matrix).  In  the  ease  of  styrene t 


Since  F is  a symmetric  matrix,  it  follov/a  that  the  eigenvalues  F are  real  (by  virtue 
of  (2.49)  they  are  « 2],  and  that  the  eigenvectors  belonging  to  different  eigenvalues 
are  mutually  orthogonal.  Prom  (2.J4)  we  know  that  the  eigenvectors  are  normalized. 
Different  eigenvectors  of  one  degenerate  level  can  of  course  be  orthogonal ized.  Hence 
the  eigenvectors  fg,  ...  of  (2.52),  corresponding  to  the  eigenfunctions  ^^(x), 
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0g(x)  ...  0^(x)  Eiatlafy  the  relations 

» N 

^np'’’mp  “ ®mn  ' 

In  analogy  with 


(2.54) 


{ " «mn  * (2.55) 

I 

The  eigenvalues  are  found  as  solutions  of  the  seoular  equation  oorrespondlng  to 
(2.52)  or  (2.47)i 


IP  - Pl|  - 0 , (2.56) 

or 

I Pi-  - 0 . (2.5'h 

Prom  the  eigenvalues  P^^  one  finds  the  energy  levels  by  means  of  (2.49)  and  (2,14'). 

When  the  eigenvalues  are  found,  the  eigenveotors  oan  be  determined  as  solutions 
of  the  system  of  equations  (2.52),  or  the  eigenvector  oan  be  determined  as  solutions 
of  (2.47).  Finally,  the  free-eleotron  wave functions  ♦^^(x)  oan  be  calculated  from  the 
eigenvectors  by  the  method  outlined  In  Eqs.  (2.21-30). 

In  this  way,  energy  levels  and  eigenfunctions  of  the  free-electron  problem  with 
constant  potential  oan  be  determined  without  explicit  use  of  the  boundary  condition 
(1.34)  and  the  Joint  oonditlons  (2. 4, 5, 6))  and  application  of  the  model  to  larger 
moleoules  Is  made  possible.  Whereas  the  determination  of  the  eigenvalues  must  always 
be  done  by  means  of  the  seoular  equation  (2.56,57)  (except  In  the  very  simplest  oases). 
It  Is  sometimes  more  convenient  to  determine  the  free-electron  wavefunotlons  dlreotly 
than  to  calculate  the  eigenveotors.  This  Is  the  case  when  the  phases  of  the  wave- 
function  are  immediately  given  by  the  symmetry  of  the  system  (for  example  In  naphtha- 
lene); once  the  eigenvalues  and  all  phases  of  a wavefunctlon  are  known.  It  Is  easy 
to  determine  its  amplitudes,  a^,  by  use  of  the  equations  (2.4)  and  (2.7). 

It  should  be  noted  that  the  matrix  form  of  the  eigenvalue  problem  oan  equally 
well  be  developed  for  eigenvectors  of  the  type  discussed  in  Eq.  (2.4o).  A nice  appli- 
cation of  this  fact  to  the  stilbene  molecule  of  Figure  7 will  be  dlsoussed  later. 
Another  consequence  la  that  the  orthonormallty  relations  (2.54)  hold  also  for  the 
eigenveotors  (2.40).  Incidentally  It  may  be  mentioned  that  the  orthogonality  of 
eigenvectors  belonging  to  different  eigenvalues  can  also  be  derived  directly  by  the 
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methofj  employed  In  the  proof  of  the  normalization  theorem  (2.34). 
d.  i piaouBBlon.  ApplioatlonB 

The  N values  which  are  obtained  for  P «•  2ooa/c  give  rise  to  Infinitely  many  values 
of  K.  There  are  N values  /e^(n  ••  1,  2...N)  In  the  range  0 < < tt  and  each  of  these  N 

values  )Cj,j  gives  rise  to  an  Infinite  set  of  k values,  namely 

*nv  “ *''n  ‘ V « 1,  2,  3 ...  . (2.58) 

A possible  situation  In  the  case  N ■ 2 is  ahovm  In  Figure  8.  We  note  that  by  defini- 
tion (c  « 0 (see  2.12).  To  each  value  = *)  belongs,  by  means  of  (2.l4'),  one 

eigenvalue  and,  by  means  of  Eqs.  (2.21-30),  one  wavefunetion  ^*^(x). 

Thus  each  eigenvector  gives  rise  to  an  Infinite  set  of  elgenfunotlons  O^^(x). 
For  the  N lowest-enargy  levels,  i.  e. , those  In  the  range  (0  « k < tt),  the  de  Broglie 
wavelength  (2. 13)  Is 

. ' (2.59) 

9 -I 

Hence  there  oan  be  no  more  than  one  node  from  each  wavefunetion  between  two  neighboring 
atoms.  For  the  wavefunotlona  corresponding  to  k > w,  on  the  other  hand,  there  always 
exists  pairs  of  neighboring  atoms  between  which  are  at  least  two  nodes. 

The  eigenvalue  ° (corresponding  to  Pj^  - 2)  Implies  a »,  i.  e. , a constant 

wavefunetion.  Hence  it  does  not  occur  if  the  electron  path  has  a free  endpoint,  since 
a constant  wavefunetion  would  have  to  vanish  Identically.  This  eigenvalue  always 
occurs,  however,  if  the  electron  path  has  no  free  endpoints. 


— 1 1 II  I U-— I 

0 27T  47T  67T  87T 

Fig.  8 

Since  the  molecule  contains  N Tr-electrons , the  one-half  N lowest  levels  are  doubly 
occupied  In  the  ground  state.  If  N is  odd,  then,  the  highest  level  Is  only  singly 
occupied  In  the  ground  state.  Later,  It  will  be  shown  that  for  alternant  hydrocarbons 
all  ground-state  levels  have 

K <s  7t/2  ; (2.60) 

but  also  for  those  non-alternant  hydrocarbons  which  were  investigated  by  one  of  us,^ 
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the  condition  (2.60)  was  always  found  to  be  valid  for  the  ground-state  levels. 

The  exoltatlona  of  the  highest  electrons  to  the  next  higher  levels  correspond  to 
the  TT-eleotron  spectra  of  the  conjugated  molecule.  The  transition  energy  AE  = Ej^^  - 
between  two  levels  E^,  E^  Is  given  by 

AE  - Ej3A((C^)  , (2.6l) 

with 


the  corresponding  wavenumber  (l/X)  •>  AE/ho  by 

(lA)  - (1/Xjj)A(k®)  , 
or 

(lA)  » iXoAdA®)  , 

with 

(1/Xjj)  « Xg/8ffV  . 


(2.62) 


(2.6J) 


(2.64) 


(2.65) 


Here  the  symbol  Ax  always  has  the  meaning  x^,  - x^^j  Is  the  de  Broglie  wavelength 
(2.13)  of  a free-eleotron 'orbital j a o 0.529151^  Is  the  Bohr  radiusj  Ejj  • 13.6035eV  Is 
the  atomic  Ionisation  potential  of  hydrogenj  and  X^  ■ 0.0242607?  is  the  Compton  wave- 
length. For  D « Dq  ■ 1.40000?  (the  subscript  C stands  for  "conjugated",  the  value 
1.40000  being  a mean  value  for  the  bond  distance  In  conjugated  systems)  one  obtains 
the  msnerloal  values 


«»  1.9432eV  a 44.827koal/mole  (2.66) 

(1A(.)  - 15,677cm"^  . (2.66') 

For  .other  values  of  D one  finds  and  X^  by  means  of  the  relation 

(Ejy/Ec)  - (X^/Xp)  - (D(j/D)2  . (2.67) 

Because  of  Eq.  (2.60)  It  follows  then  that  the  lowest  7r-eleotronlc  excitation  energies 
are  of  the  order  of  magnitude  given  In  Eq.  (2.66,66'),  or  rather  a little  larger.  This 
result  Is  In  agreement  with  experiment.  In  contrast,  the  excitation  o.«'  an  electron 
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p Q 

from  a ground-state  level  to  a level  with  /c  ? tt  needs  an  energy  AE  » (rr  - ” 

l4.j5eV.  Since  the  ionization  energy  for  benzene  Is  9«2^eV,  such  a level  cannot  have  a 
physical  meaning,  1.  e.,  the  approximations  of  the  free-electron  model  which  do  not 
allow  for  ionization  become  insufficient  for  excitations  of  this  magnitude. 

The  lowest  N levels  of  the  free-electron  model  correspond  bo  the  N levels  found 
in  LOAD  theory.  The  similarity  between  the  form  of  the  matrix  F (2.53)  and  the  matrix 
which  determines  the  eigenstates  In  LCAO  theory  la  strlklngj  indeed,  the  only  differ- 
ence Is  that  each  joint  Introduces  a factor  (2/3)^  Into  the  matrix  F,  whereas  It  does 
not  do  so  In  the  LCAO  matrix.  The  set  of  linear  equations  (2.44,45,46)  Is  of  the  same 
nature  as  the  corresponding  system  In  LCAO  theory^  and  the  slmpllfluatlons'  on  the  basis 
of  the  molecular  symmetry  can  therefore  be  achieved  by  exactly  the  same  methods  which 
are  used  in  LCAO  theory.^®  It  should  be  noted  that  each  free-electron  eigenfunction 
belongs  to  a definite  representation  of  the  three-dimensional  symmetry  group  of  its 
molecule,  In  order  to  find  this  representation  one  has  to  take  Into  account  that  the 
total  three-dimensional  free-electron  wavefunctlon  (see  Eq.  (1.17)  ff.)  IB  anti- 
symmetric with  respect  to  the  molecular  plane.  It  follows  that  a one-to-one  corres- 
pondence exists  between  the  free-electron  states  and  the  LCAO  states  according  to  their 
group-theory  classification.  However,  the  order  of  the  energies  Is  not  alwSys  the  same 
same;  an  example  is  given  by  the  fourth  and  fifth  levels  in  perylene  (of.  refs.  6 and 

8) I The  LCAO  model  furnishes  E.  > E-  , the  free-electron  theory,  Ep  > E,  , 

*u  ^lu  “lu  *u 

Also  the  free-electron  method  may  yield  accidental  degeneracies  not  shown  by  the  LCAO 

model  and  vice  versa.  (E.  £. , the  sixth  to  ninth  levels  inclusive  In  perylene  are 

degenerate  In  the  LCAO  model,  but  not  in  the  free-electron  model.) 

Let  M_  and  be  the  two  bond  midpoints  on  both  sides  of  a non- Joint  atom  P,  and 

let  P and  Q be  the  two  atoms  on  both  sides  of  a bond  midpoint  M.  Then  it  is  sensible 

to  define  an  atom  population,  a',  on  P and  a bond  population,  b' , on  the  bond  PQ  by 

the  equations 

a'(P)  - E fj^  ; dx0^(x)  , (2,68) 

n M_ 

b»(M)  - b'(PQ)  " J fn  I ' (2.68') 

where  fj^  ■ 0>  2 denotes  the  number  of  electrons  occupying  the  orbital  0^(x).  The 

sum  extends  only  over  the  filled  orbitals,  since  f ^ a 0 for  the  other  cases.  For  an 

^®See,  e.  £.,  ref.  8,  page  I96. 
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atom  Pj  at  a Joint  the  atom  population  has  to  be  defined  by 

3 Mr 

J ^n  gf ^ f (2.68") 

J 

where  Mg,  Mj  are  the  three  bond  midpoints  between  Pj  and  ite  three  neighbor  atoms. 

Instead  of  using  the  quantities  a'  and  b«.  It  Is  simpler  to  work  with  the  quan- 
tities a and  b,  defined  by  (See  Eq.  (2.32)) 

a(P)  - L fjjbnp  - m(P)  , (2.69) 

n 

b(Pft)  - b(M)  - D tu(M)  , (2.69' ) 

where 

o)(x)  - Sfn0n(x)  (2.70) 

denotes  the  total  m-eleotronlo  linear  density  at  the  point  x,  and  Tp  Is  defined  in 
Eq,  (2.50).  One  would  have  a ■ a'  and  b - b',  If  ♦(P),  ^(M),  ^(Pj)  would  represent 
the  mean  value  of  the  function  d(x)  in  the  integrals  (2.68),  (2. 68'),  (2.68"),  respeo- 
tlvely.  Although  this  Is,  in  general,  not  true  for  a single  atom  or  bond.  It  oan  be 
considered  to  hold  ^ the  average,  in  the  sense  of  the  following  argument.  According 
to  Eqs.  (2.33)  and  (2.43)  the  total  TT-electronie  population  of  the  molecule  oan  be 
expressed  as 


whence 


N - E f / dx*2(jt)  » £ a(p)  - S b(M)  , 
n " J " P M 


(2.71) 


E a'(P)  o r a(P)  « N , (2.7I') 

P P 

E b'(M)  » E b(M)  = N , (2.71") 

M H 

where  the  summation  over  P oontalns  all  atoms  and  the  summation  over  M contains  all 
bond  midpoints  (N  Is  the  number  of  electrons).  Hence,  on  the  average.  a(P)  may  be 
considered  as  the  contribution  of  atom  P to  the  total  population,  and  b(PQ)  may  be 
considered  as  the  contribution  of  the  bond  PQ  to  the  total  population.  In  the  follow- 
ing, we  shall  apply  the  terms  atom  population  and  bond  population  to  a and  b respec- 
tively. 

The  quantities  a,  b have  a similar  physical  meaning  as  the  analogous  quantities 
q («  the  number  of  electrons  on  an  atom)  and  p (■  bond  order  between  two  atoms)  of  the 
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LCAO  model. 2 it  follows  then  that  the  reactivities  of  the  different  atoms  and  the 
total  7T-eleotronlo  dipole  moment  of  the  molecule  are  connected  with  the  a 'a  In  the  same 
way  as  they  are  connected  with  the  q's,  and  the  application  can  be  taken  over  directly 
from  LCAO  theory,  it  should  be  noted  however  that  at  the  Joints  there  is  a difference 
between  the  electron  density  and  the  electron  population  (f  Is  difference  Is  unknown  In 
LCAO  theory),  and  one  may  i’alse  the  question  whether  the  former  or  the  latter  decides 
the  chemical  reactivity.  OTie  quantities  b,  on  the  other  hand,  are  related  to  the  bond 
lengths  and  the  so-called  free  valences  of  the  same  atoms  In  the  same  way  as  the  LCAO 
quantities,  p. 

On  this  basis  It  is  possible  to  correlate  the  free-electron  wavefunotlons  with 
experimental  facts,  The  numerical  comparison  shows  that  the  free-electron  results  are 
In  no  way  inferior  to  the  results  obtained  from  the  LCAO  model, 
e.  Alternant  Conjugated  Systems 

We  conclude  the  present  paper  by  showing  that  the  free-electron  theory  describes 

the  particular  character  of  the  alternant  conjugated  system  with  an  even  number  of 

11 

atoms  in  the  very  same  way  as  the  LCAO  model  does.  This  important  group  of  mole- 
cules, which  includes  all  conjugated  structures  composed  of  even-membered  rings,  Is 
defined  by  the  following  properties i The  conjugated  structure  contains  an  even  number 
of  atoms , and  ^ must  be  possible  ^ divide  all  atoms  Into  two  classes . e , £. , "starred 
ones"  and  "uns tarred  ones".  In  such  a fashion  that  no  two  neighbors  belong  to  the  same 
class. 

The  theoretical  characteristics  of  these  molecules  have  their  origin  in  the 
"pairing  theorem",  which  we  shall  now  establish.  Let  J be  the  matrix 


r 1 0 0 0 ...  0 1 

0-1  0 0 ...  0 


L 0 0 0 0 ...±i  J 

1.  e.,  J 19  a diagonal  matrix  In  which  all  elements  corresponding  to  the  starred  atoms 
are  1 and  all  elements  corresponding  to  the  unstarred  atoms  are  -1.  Then 

= I , (2.73) 

^^Alternant  systems  with  an  odd  number  of  atoms  present  no  difficulties  and  v;lll  be 
dealt  with  In  a subsequent  paper. 
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where  I denotes  the  unit  {natrlx,  tind 

JFJ  = -P  , 

where  P Is  the  matrix  (2.53).  Henoe  It  follows  from  (2.52)  that 

p^i  ™ pi^i  , 

with 

pi  - -p  , 


(2.74) 

(2.75) 

(2.76) 

(2.77) 


That  Is:  If  ^ Is  an  elgenveotor  with  the  eigenvalue  F,  then  Is  an  eigenvector  with 

n 

(2.78) 


the  eigenvalue  F>  <•  -F.  Henoe,  If  Is  an  eigenvalue  of  k (see  Eq.  (2.49)),  then 


n - 


Is  also  an  eigenvalue.  Combined  with  Eq.  (2.58),  this  leads  to  the  result  that  with 
/c^(<  7t/2)  one  has  simultaneously  the  eigenvalues 


«nv  - + vn 


*nv  “ 


1 V = 0,  1,  2,  3 
7t-K«)  + V7r  J 


(2.79) 


Figure  9 shows  all  the  eigenvalues  generated  by  one  k^,  which  must  appear  together  in 
an  alternant  system. 


— 1 

1 ■■ 



—4- 

> » 

27T 

FlS.  9 


37T 


41T 


Let  be  the  smaller  one  and  let  be  the  larger  one  of  the  two  paired  eigen- 


values « 7T,  1..  e . , 


0 < « #7T  , 


^7T  € « 7T  . 


(2.80) 

(2.81) 


Then  it  follows  for  the  two  paired  eigenfunctions  0^(x)  and  ^,!j(x)  that 
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, (S.82) 

4D  ?>  & 2D  , (2.83) 

i.  e. , for  ♦n(x)  two  nodes  are  at  least  two  bond  lenistha  apart,  and  for  d^(x)  the 
distance  of  two  nodes  lies  between  one  and  two  bond  lengths.  Moreover,  one  sees  from 
Eq.  (2.76)  that 

starred  atom,  . 

I (a. 84) 

^n^^^  ■ at  each  unstarr.ed  atom,  J 

or  vlee  versa,  depending  on  the  oholoe  of  the  arbitrary  signs  in  the  wave funot ions. 

The  explicit  form  of  the  wavefunotlon  ^^(x)  becomes  simple  on  a branch  where  the 
point  Xg  - 0 oolnoldos  with  an  atomi  If  ^^^(x)  given  by  Eq.  (2.2),  then  ^^(x)  is 
given  by 

“ ®Bn  ao8[(D‘Vk„)Xs  - 63^]  , (2.85) 

as  may  be  verified  by  taking  account  of  Eqs.  (2.78,84)  in  Eqs.  (2.23,  24)  which  deter- 
mine kj^,  k^,  6^^.  From  (2.85)  follows  indeed  that 

“ (-l)''*Bn(''°)*  V - 1,  2,  3 , (2.86) 

which  is  identical  with  (2.84),  if  one  chooses  the  point  « 0 to  be  a starred  atom. 

On  the  other  hand,  on  a branch  where  Xg  * 0 coincides  with  a bond  midpoint  one  finds, 
by  virtue  of  (2.29,30)  and  (2.78,84),  that 

" ®Bn  oos[(D"^-k„)xg  + (^TT-6g^)]  , (2.87) 

whence 

«tin[(v-i)D]  - (-!)%„[  (v-^)D),  V = 1,  2,  3 ...  , (2.88) 

which  agrees  with  (2.84)  if  Xg  a -^D  Is  chosen  as  a starred  atom. 

An  immediate  consequence  of  the  pairing  theorem  Is  the  population  theorem,  which 
states  that  the  total  7r-electronic  population  at  each  atom  is  unity  in  the  ground 
state  of  an  alternant  conjugated  system.  This  theorem,  which  excludes,  for  example, 
the  existence  of  a 7r-eleotronlo  dipole  moment  for  alterr.ar.t  systems,  is  established  in 
the  following 'way : By  virtue  of  the  pairing  theorem,  in  particular  Eq.  (2.84),  the 


I, 

t 
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TT-eleotronlo  atom  population  (2.70)  in  the  ground  state  can  he  written 

N/2  5 N « 

a(P)  » r 2<p2p  - S (p2p  . (2.89) 

noi  n=l 

Now  Eq.  (2.64)  states  that  the  quantities  f^^p,  taken  for  n o i,  2,  . . . N,  form  an 
orthogonal  matrix  and  hence 

a(P)  - 1 . (2.90) 


This  population  theorem  Is  the  counterpart  of  the  Identical  theorem  In  the  LCAO  model. 

Since  the  secular  equation  slmultaneovualy  has  the  eigenvalues  P and  -F«  It  must 
be  possible  to  write  It  as  a system  of  equations  containing  only  P . Such  a system 
can  be  obtained  in  the  following  mannert  The  matrix  ooiranutss  with  J (2.72) s |i; 


- J(-r)(-D  - , 


(2.91) 


and  hence  has  no  elements  (different  from  cero)  connecting  starred  and  unstarred  atoms. 
Hence  the  vector  transformation  described  by  F leaves  both  the  "starred  subspaoe"  and 
the  "unstarred  subspaoe"  Invariant,  and  hence  the  eigenvectors  of  T o&n  be  assumed  to 
have  only  starred  components  or  only  unatarred  components.  The  eigenvalue  of  is 
of  course  P^.  Indeed  from  (2.52)  and  (2,75»  76>  77)  it  follows  that 


y®#  - p®#  » 

FJ#  - P®  Jfl  , 

whence 

with 


= id  i j)# . 


^ 1 
1 

0 

1 

1 

0 

0 

, (1  - J)  - 

1 

• 

. • . 

e 

, 

(2.92) 

(2.92») 

(2.95) 

(2.95*) 

(2.94) 

(2.95) 


I 


..'4- 
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one  Tnay  write  Eq.  (a. 95)  and  (2.95')  ae 

(F®)*  K - . (2-96) 

(^^0  *0  " , (2.96') 

whore  la  the  4N-dlmenalonal  veotor  obtained  from  by  omitting  all  unetarred 
components,  and  is  the  iN-dlmenslonal  veotor  obtained  fr-om  by  omitting  all 
starred  components.  Similarly,  (7®),,  and  (7®)^  are  obtained  from  7®  by  omittli:®  all 
unstarrert  or  all  starred  elements  respectively.  The  eigenvalues  P®  are  therefore 
solutions  of  the  two  equations 

1(7®)*  - P®1|  » 0 , (8.97) 

|(7®)o  - P®1|  - 0 . (2.97') 

The  eigenvectors  of  (2.96)  and  (2.96'),  corresponding  to  the  same  eigenvalue  F®, 

furnish  the  two  elgenveotors  and  of  7®.  Prom  these,  one  obtains  the  two  eigen- 
vectors 

(2.98) 

#'•  . (2,98') 

of  7i  corresponding  to  the  two  eigenvalues  +T  and  -F.  VThloh  eigenvalue  belongs  to 
which  eigenvector  is  decided  by  checking  the  relation  (2.44),  1.  e. , one  of  the 
equations  of  the  set  (2.52),  between  three  neighboring  points. 

As  an  Illustration  let  us  finally  consider  the  problem  of  finding  the  eigenvalues 
of  the  stllbene  molecule,  whose  free-electron  path  has  already  been  shown  in  Figure  7. 
In  the  first  place,  we  note  that  one  may  set  up  a secular  equation  for  the  four  points 
A,  B,  C,  D,  which  liave  the  neighbor  distance  5D.  Since  they  fom  an  alternant  system, 

A* 

the  four  eigenvalues,  P^,  of  this  problem  form  two  pairs,  and  since  there  are  no  free 
endpoints,  one  pair  must  be  P^^  « 2,  P[^  « -2,  corresponding  to  = 0,  = it,  where 

Ic  a 5Dk.  Hence  there  actually  remains  one  eigenvalue  to  be  determined.  The  matrix 
(2.55)  becomes  for  these  four  points. 
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B 


Whence 


A 

B 

C 

D 


2(f) 

0 


8 

5 


TV 


2{|)^ 

0 

2 

7 

0 


2 

5 

0 

2,i 


2(f) 


28 

T 


^(S\i 

TV 


^(|)i 

0 


0 

0 


2(f) 


4,2vi 

TV 


8 

? 


and 


(P^)* 


28 

T 

5(5)  • 


'♦fsa 

8 


(f®). 


The  eeoular  equation  (2.97)  becomes 

(JF)*^  " 52(3P)®  + 576  o 0 , 


(2.99) 


(2.100) 


(2.101) 


(2.102) 


and,  since  (3P)®  = 36  (corresponding  to  P « i 2)  must  be  one  root,  It  follows  that 
(3P)®“  52  - 36  ■ 16  is  the  second  root.  Thus,  we  find  for  iP  ■ oosK,  and  for  k,  the 
four  eigenvalues  ' - 


oosK  I 1 2/3 -2/3  -1 

Kjj  I 0 0.2687r  0.7327r  tt 

in  the  range  {0  ^ k ^ tr) . These  eigenvalues  plus  those  generated  in  the  range  between 
7T  and  3'»f  are  drawn  In  Figure  10  on  horizontal  lines.  The  first  vertical  scale  shows 
the  values  of  k and  the  second  (middle)  vertical  scale  shows  the  corresponding  values 
of  k B k/3D. 

Since  the  eigenvalues  of  k depend  only  on  the  free-eleotron  path  and  not  on  the 
manner  in  which  one  assumes  equidistant  points  on  the  path,  no  other  eigenvalues,  k, 
can  be  found  by  solving  the  l4-dlmenslonal  secular  equation  corresponding  to  the  14 


52 


FREE-ELECTRON  MODEL.  1 


atomB  of  the  itioleouie.  However^  In  this  problem  one  works  with  the  quantities  k - Dk  ■ 
k/3>  The  14  states  with  0 k « tt  oorrespond  to  the  14  LCAO  states  of  the  moleoules. 
Kow  the  third  (right)  soale  in  Figure  10  shows  the  values  of  ye  for  the  . eigenvalues 
whloh  were  found)  one  recognises  that  there  exists  only  10  levels  in  this  range;  hence 
some  of  these  must  be  degenerate.  Indeed,  the  levels  o tt/j  and  le^  = 27t/3  are  both 
triply  degenerate,  thus  providing  the  extra  four  levels.  This  degeneracy  IS  easily 
recognized:  The  eigenfunctions  Ibr  b tr/3  have  » 6D,  and  one  sees  Immediately 
that  one  can  construct  three  mutually  orthogonal  wavefunotlons  with  this  wavelength; 
they  are  shown  In  Figure  11  where  a 

full  line  Is  drawn,  wherever  the  wave  funotion  is  positive, 
hollow  line  Is  drawn,  wherever  the  wave  fimotion  is  negative, 
dotted  line  is  drawn,  wherever  the  wave  funotion  is  Identically  zero. 

Since  there  are  three  Independent  wavefunotlons  for  k = tt/35,  the  same  must  be  true  for 
X o 27t/3  by  virtue  of  the  pairing  theorem.  In  order  to  find  all  wavefunotlons  one 
solves  the  system  of  two  linear  equations  obtained  from  (2.96)  by  substituting  (F^),j 
of  (2,101)  for  (F^)#;  by  the  method  described  after  (2.97),  one  then  obtains  the  eigen- 
vectors of  F (2.99)  and  hence  by  the  usual  methods  [see  Eqs.  (2.21-30)]  the  wave- 
funotlons for  the  14  lowest  states.  Finally  one  can  determine  the  symmetries  of  the 
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Fig.  11 


different  wavefunotlons. 

Of  course  it  would  have  been  possible  from  the  beginning  to  separate  out  the 
eigenvalue  equations  corresponding  to  the  different  symmetry  types j this  procedure 
would  be  advantageous  In  dealing  with  a molecule  of  sufficiently  high  symmetry.  Since 
this  method  operates  In  prsbieely  the  same  way  as  In  the  LCAO  model, we  do  not  have 
to  go  Into  details  at  this  point,  and  refer  the  reader  to  the  following  paper^  for  an 
Illustrative  example. 


APPENDIX 

JUSTIFICATION  OP  THE  BOUNDARY  CONDITION  FOR  FREE  ENDPOINTS 

The  assumption  that  the  wavefunotlon  0(x)  goes  to  zero  one  bond  length  D beyond 
the  last  atom  at  a free  end  [see  comments  after  Eq.  (l.;?4)]  has  been  essential  In 
deriving  the  normalization  of  the  elgenveotors  [see  after  Eq.(2.?8)]  and  In  deriving 
the  matrix  formulation  of  the  eigenvalue  problem  [see  Eq.(2.45)].  Hence,  also  the 
orthogonality  of  the  elgenveotors  depends  on  the  present  oholoe  of  the  boundary  oondl- 
tlons  [see  comments  after  (2.5?' )]• 

It  will  now  be  shown  that  this  boundary  condition  is  the  only  one  which  guarantees 
the  normalization  of  the  eigenvectors  as  well  as  the  existence  of  the  pairing  theorem 
for  alternant  systems . 
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Consider  tho  proof  starting  after  Eq.  (2.^4).  Let  us  assume  that  the  branch  B 
oonaldered  there  has  a free  endpoint j let  the  "first  atom  Pg"  [definition:  between  Eqs. 
Eqs.  (2.34)  and  (2,35)]  the  one  nearest  to  the  free  endpoint}  and  let  the  endpoint 
of  the  electron  path  be  located  at  a distance  (D  beyond  the  atom  Fg.  Then  the  branch 
length  becomes 

ig  o D(n  - 1 + 5)  . (A.l) 

We  ohoose  the  path  endpoint  to  be  the  origin  Xg  - 0 of  the  coordinate  Xg,  which 
increases. in  the  direction  away  from  the  endpoint.  Hence,  the  wavefunotlon  vanishes 
for  Xg  ■ 0,  and  has  therefore  the  form 

♦g(xg)  - ag00B(kXg  - itr)  , (A. 2) 

so  that  the  argument  of  the  cosine  becomes,  for  the  first  atom  Fg, 

f - e»c  - it  , (A.3) 

[see  Bq.  (2.35)]. 

The  course  of  the  proof  subsequent  to  Eq«  (2.37)  shows  that  under  the  present 
curcuma tances  [note  Eq.  (A.l)]  the  Eqs.  (2.33)  result  from  Eq.  (2.37)  only  If  the 
relation 

-Cota  slnf  oosf  + oos®f  ■ i (A. 4) 

is  fulfilled . Hence , 

- - 1)  • (*-5) 

and  by  virtue  of  (A.3)« 


slnqK  o n aina  , 


(A.6) 


with 


TJ  «.  24  - 1 . 


12 

If  there  are  several  branches  with  free  endpoints,  the  Eq.  (A. 4)  must 
for  each  branch  separately  In  order  that  our  results  remain  generally 


(A. 7) 


still  hold 
valid. 
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Eq.  (A. 6)  Is  satisfied  by  the  values 

Tj  ■ 0 I +1  » <■!  » (At8 

Now  there  are  always  eigenfunctions  with  k In  the  range  0<ic<frr  (k«0  Is  excluded 
If  there  is  a free  endpoint).  For  such  values  of  k there  exist  no  other  solutions  of 
(A. 6),  as  may  be  seen  In  the  following  wayt  Figure  12  shows  the  two  curves 


f(n)  B RlnuK  , 


(A. 9) 


s(n)  **  n Blnx  . 


(A. 10) 


Fig.  12 

For  Tj  o (tt/2>c)  the  straight  line  (A. 10)  assumes  the  value 

g(n/2K)  =»  (it  8lnK)/2>c  , 


which  decreases  from 


'V 


JlSl.  - 


; - ^v■.. 


Hr..,. 
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g(7r/2K)  * 1 , for  K ■ Jit  . 

Hence  the  line  (A. 10)  must  lie  in  the  sector  between  the  two  dotted  lines  of  Figure  12; 
as  an  example  the  curve  for  k « n/3  has  been  indicated.  (It  should  be  noted  that,  as 
K -•  0,  the  argument  point  tt/Sk  moves  out  towards  Infinity.)  It  follows  then  that  the 
two  curves  (A. 9)  and  (A. 10)  have  no  other  intersections  than  the  ones  given  by  (A. 8). 
Thus,  Eq.  (A. 5)  is  only  fulfilled  by  (See(A.6)] 

e - 0 , J , 1 . (A. 12) 

Since  the  wavefimctlon  shall  certainly  not  vanish  at  the  last  atom,  the  first  value 
t a ;o  has  to  be  discarded.  The  second  value  ( a J is  Inoonslstent  with  the  validity 
of  the  pairing  theorem.  Indeed,  it  la  easily  recognised  that  such  a boundary  condi- 
tion would  lead  to  the  relation 

- (1  + 2oo8/e)«3(PQ)  - 0 (A. 13) 

Instead  of  to  Eq.  (2.45).  Consequently,  the  matrix  f entering  in  the  eigenvalue 
equation  (2.51)  would  contain  the  diagonal  element  (-1)  for  each  end  atom;  and  hence 
Bq.  (2.74),  the  basis  for  the  pairing  theorem,  would  not  hold. 

Therefore  the  normalization  of  the  eigenvectors  and  the  pairing  theorem  for 
■ alternant  systems  result  only  from  the  choice  $ ■ 1;  and  this  is  the  definition  of  the 
path  end  point  adopted  in  this  paper. 
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